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SUMMARY 

We simulated the remagnetization dynamics of the ultra-dense and ultra-thin magnetic dot array system with 
dipole-dipole and exchange coupling interactions. Within the proposed 2D XY superlattice model, the square dots 
are modeled by the spatially modulated exchange-couplings. The dipole-dipole interactions were approximated by 
the hierarchical sums and dynamics was reduced to damping term of the Landau- Lifshitz- Gilbert equation. The 
simulation of 40000 spin system leads to nonequilibrium nonuniform configurations with soliton-antisoliton pairs 
detected at intra-dot and inter-dot scales. The classification of intra-dot magnetic configurations was performed 
using the self-adaptive neural networks with varying number of neurons. 
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1. INTRODUCTION 

Many interesting applications of the periodic 
magnetic dot arrays are expected in the field 
of the magnetic recording [1] and magnetic sen- 
sors. The typical properties brought by the ar- 
ray geometry is the limited length of the spin 
chains coupled by the exchange coupling and the 
uniqueness of the magnetostatic interactions at 
the scales comparable or larger than dot size. 

The powerful tools to study the static and dy- 
namic properties of magnetic nanostructures in- 
cluding the thin films and isolated small particles 
represent micromagnetic [2], and Monte Carlo 
[3, 4, 5, 6] simulations. They are focussed to the 
analysis of the nonuniform magnetization states. 
On a purely qualitative level the intra-dot nonuni- 
formities can be characterized as a mixtures of 
vortices, flowers, domain walls, etc.. Our recent 
activity is related to the search for the tools allow- 
ing the systematic classification of the simulated 
magnetic structures. 

The methodology presented in this paper is 
inspired by the progress in the theory of artifi- 
cial neural networks, that are nonlinear models 
suitable to reduce, classify or interpolate the in- 
formation involved in the systems of the complex 
patterns [7]. In the previous paper [8] we sug- 
gested the implementation of radial basis function 
networks to model the magnetic dot array mag- 
netization processes. The lack of this proposition 
was the absence of the network learning which 
needs the data support from physically accept- 
able micromodel simulations. 

The above reasons led us to make some 
progress towards the opposite direction by us- 
ing the neural network analysis of the configura- 
tions generated by the remagnetization processes 
described by the XY superlattice model of the 
ultra-dense array of the square dots limited to in- 



plane spin polarization. The motivation for the 
study of ultra-dense spin structures is the special 
interest for the interfacial inter-dot effects. The 
experimental studies of the magnetization rever- 
sal are generally not in agreement with the fa- 
miliar model of coherent rotation. Instead, com- 
plex processes such as curling and buckling are 
observed. The model structures were obtained 
by simulations based on the numerical solving 
of Landau-Lifshitz-Gilbert equation. The simu- 
lations bring the usual problem of effective data 
post-processing [9] . In the present paper the data 
were treated by means of so called adaptive re- 
sonance neural network ART, which is introduced 
in section 3. This tool was used to classify the 
intra-dot magnetic configurations, which are the 
transients of the magnetization reversal process. 
The ART choice is connected with its ability to 
vary the number of neurons (network topology) 
according to complexity of treated inputs, which 
are in our case activated by the intra-dot mag- 
netic configurations. In the section 4 we concen- 
trate to the analysis of data provided by the simu- 
lated magnetization reversal. 

2. XY MODEL OF ULTRA-DENSE 
ARRAY 

The details of the remagnetization dynamics 
are studied by means of the model of finite mag- 
netic superlattice system described by the XY 
spin Hamiltonian 
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The first sum of Eq.(l) represents the dipolar con- 
tribution controlled by the constant K, where 
S« (|S Q | = 1) is the planar spin variable and 
is the demagnetization tensor of superlattice; 
a, $ are site indices. The second sum represents 
the exchange energy contribution with periodi- 
cally modulated nearest-neighbor exchange cou- 
pling interaction Je Q ^, where e Q/ 3 G {0, 1} (see 
Fig. 1). The third sum represents standard Zee- 
man term where e x is the Cartesian unit vector 
oriented along the main axis of superlattice. Pre- 
sented model belongs to the idealized structure, 
where the patterning attains the maximum geo- 
metric filling independent of dot size (L linear 
size of dot) . 




Figure 1: The schematic plot of the exchange cou- 
pling (bonds represents e a p = 1) modulation. The 
system contains 16 dots, each of them consists of 25 
spins. 

The model dynamics was described by the 
Landau-Lifshitz-Gilbert equation, that in 2D case 
reduces to relaxation term. The numerical in- 
tegration scheme was derived from the recursive 
propagator form 



S a (t + At) = U 
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where oj a (t) = (S a (t) ■ e x )(h a (t) ■ e y ) - (S a (t) ■ 
e y )(h a (t) -e x ) is the instant angular frequency of 
the spin rotation, h a (t) = —5H/5S a is the local 
effective field; t is the time, At is the integration 
step. The unitarity of the spin rotation matrix 
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implies that numerical integration scheme con- 
serves the spin vector size independently of ad- 
ditional approximation of integral in Eq.(2). The 
computationally demanding task of the integra- 
tion is the calculation of the dipolar fields, where 
essential reduction brings the hierarchical sum- 
mation [10]. In the case of array, the most natu- 
ral hierarchy level choice is the block association 
with the single intra-dot moment. 



3. THE ART CLASSIFICATION OF 
INTRA-DOT CONFIGURATIONS 

The artificial neural network models are in- 
spired by the physiology, and mimic the neurons 
and synaptic connections of a brain. They can 
be considered as mappings [7] constructed from 
given activation functions. The unknown param- 
eters of these functions called synaptic weights 
are adjusted by training. The fascinating feature 
of the neural systems is the associative recog- 
nition of complex structures. In this paper we 
deal with the ART networks developed by Car- 
penter and Grossberg [11], originally as a model 
to explain the adaptive phenomena in visual sys- 
tems. The family of ART algorithms belongs to 
group of unsupervised learning algorithms based 
on the theory of adaptive resonance motivated by 
the need to construct network sufficiently flexible 
to novel inputs and preserving previously learned 
patterns. ART model represents simple clus- 
tering algorithm, which has been complemented 
with the ability to generate new neurons if neces- 
sary. This is done by using a so called vigilance 
parameter. For the summary of the key results 
and examples of ART network applications see 
e.g.[12] 

These properties led us to the opinion that 
ART network should be powerful to classify and 
compress the magnetic intra-dot configurations 
generated during the reversal. First, we must en- 
code the magnetic structure into network input 
format. We used encoding, where each dot mag- 
netization field is replaced by N c effective mag- 
netic moments located at interaction centers. The 
moments income into 2N C — dimensional vectors 
[8] named as frij. where i is the dot label. The 
vectors m, represent inputs of ART which com- 
presses this information into iV w output neurons 
named uij. Their structure is 

rhi = [mji,mj2, • • • ,m iJVc ] 

Wj = [Wji,Wj2,...,Wjjv c ] 

where i G {1, 2, . . . A d }, j G {1, 2, . . . N w } and 
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index n enumerates the spin blocks; m in repre- 
sents locally averaged dot microstate over the N$ 
spins (see Fig. 2) belonging to the square element 
□i„ of ith dot. The match between ART input 
rhi an d output Wj is measured by the " magnetic" 
Euclidean distance 
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(T is superscript denotes the transposition). 
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Figure 2: The example of the averaging from Eq.(4) 
showing the sixth square segment n i6 of ith dot in- 
cluding iVs = 49 spins. 

The learning algorithm is described in the follow- 
ing points 

1. Initialization setting N w = 1, u>i = rhi for 
random selection i S {1,2, ... N^}. 

2. Loop through the set of the intra-dot mag- 
netic patterns. For each randomly selected 
dot i follows: 

2.1 Presentation rhi to the network. The learn- 
ing is repeated for iteration steps k = 
1,2,... 

2.2 Computing of the actual index j*(k,i) of 
the winning neuron Wj*(k,i)(k) trained by 
rhi according to the competitive rule 

j*(k,i) = arg min \\uij(k) — rhi\\ 

i=l,2,...,JV w 

(6) 

2.3 Comparing \\wj*( ki )(k) — rhi\\ to the vigi- 
lance parameter p. 

2.3.1 Update of the weights via the Hcbbian learn- 
ing rule [7] 

w r{k S)(k + 1) = w r(kyl) (k) + 

rj(k) (rh t - w r{k ^(k)) , (7) 

is applied if — rhi\\ < p with the k- 

dependent learning rate parameter rj(k) = 

?70 CXp (-fc/riearn)- 

2.3.2 Creation of the new neuron wjv w +i <— rhi, 
A w <- N w + 1 if \\w r{k ^(k) -rhi\\ > p. 

2.4 Annihilation of the neuron pair 
{w z *(k),w z »(k)], z\ < z% selected ac- 
cording to relation 

[z*,Z2\= arg min \\w Zl - w Z2 \\ 

Zl,Z 2 

(8) 

if \\w z * - w z * || < p. 

The product of annihilation [z^zj] — > z\, 
N w <— N w — lis the neuron w z * determined 
by the midpoint rule w z * <— i (w z * + w z * ) . 
After w z * update Wj> z * neurons undergo to 
collapse: wj-\ <— Wj. 



3. Stop criterion is represented by the inequal- 
ity 

l N d 

(9) 

where e is small parameter. If the above 
inequality is not fulfilled the algorithm fol- 
lows from the step 2 with k incremented by 
1. 

The fulfillment of the last criterion means that 
network attains a fixed point. By this way fixed 
neurons uij, j = 1,2,..., A^ w represent the collec- 
tion of the basic types of the extracted intra-dot 
configurations. In the case of magnetization re- 
versal, the ART network has been applied inde- 
pendently to arrays at different time steps. The 
typical intra-dot and inter-dot remagnetization 
transient states are presented in the next section. 

4. RESULTS OF SIMULATION 

The computations were performed for system 
of 100 dots which each dot consists of 400 spins 
in external field H/J = — 1 and K/J = 0.1 with 
initial saturated state M(t = 0) • e x = 1, where 
M is the magnetization of array. 

In the array simulations we observed the ho- 
mogeneously magnetized dots as well as dots 
showing the nonuniform remagnetization. At the 
beginning of reversal the opposite external mag- 
netic field forms domain nucleation centers at the 
dot corners. The domains inclined by the exter- 
nal field then grow and pass towards the dot cen- 
ters by diminishing the area of the opposite do- 
mains. If the same polarity domains originating 
at different parts (corners or edges) of the dot are 
separated by the domain of opposite polarity, the 
domain growth causes the formation of two-fold 
domain wall - so called soliton-antisoliton pairs 
[3]. 

The homogenization of intrinsic dot region is 
attained by the driving effect of external mag- 
netic field which leads to wall motion annihilating 
the soliton-antisoliton pairs (see model Fig. 3). 
The coherent rotation occurs only at the early 
stages of the remagnetization for the sufficiently 
small dots when the exchange length is suffi- 
ciently greater than size of dot. This coherence 
locates nearly the corners of array. The formation 
of the soliton-antisoliton domain walls is typical 
for the central parts of the array and later stages 
of the remagnetization. 
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Figure 3: Two schemes of the principal tran- 
sient states of the nucleation, soliton-antisoliton pair 
formation, domain motion and double-wall annihi- 
lation. The axial-symmetric scenario with open- 
wall soliton-antisoliton pair (case a) is typical for 
the earlier transients of remagnetization and non- 
centric dots, whereas the center-symmetric soliton- 
antisoliton closed- wall mode (case b) is typical at the 
later transient states and dots located near to the 
center of array. 

The main point of analysis of intra-dot magne- 
tization configurations deals with the implemen- 
tation of ART network treating the subsequent 
stages of the remagnetization. The ART training 
was realized for the configurations at times sepa- 
rated by the time step At for training parameters 
770 = 0.1, T lcarn = 20 and e = 10~ 6 , p = 1.15. The 
optimized choice of vigilance (p = 1.15) stems 
from the preliminary quasistatic calibration simu- 
lations (see Fig. 4), where vigilance was slowly de- 
creased from the initial large value p = 1.5 (stabi- 
lizing one neuron) to the small values pushing the 
network to the incorporation of many neurons. 
The view point mediated by ART leads to the 
intra-dot taxonomy depicted by Fig. 5, where the 
nearest neurons [in the sense of distance Eq.(5)] 
belonging to different times arc connected by ar- 
rows. The ART model clearly distinguishes be- 
tween coherently rotating monodomain dots and 
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Figure 4: The pruning of ART network with the 
increasing vigilance parameter for K/J = 0.1 and 
I/O = 20. The intra-dot configuration is locally aver- 
aged over Ns = 4 spins [see Eq.(4)]. a) the configu- 
ration near to the switching time (t ~ 850 At); b) for 
the time when internal energy attains its maximum 
(t ~ 1100 At). 

dots including complex domains structures. The 
increasing diversity of configurations occurring 
immediately after the switching time is self- 
adaptively reflected by the extended population 
of neurons. The ART compression for later re- 
versal transient states is clearly caused by the 
small intra-dot diversity within the array. The 
details can be studied by means of the parame- 
ter p, which represents the virtual "magnifying 
glass" capable to focus to the most interesting 
final intra-dot transients (see Fig. 6 for p = 0.8). 
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Figure 5: The ART network representation of the magnetization reversal for K/J = 0.1, Lo = 20. The con- 
figuration belongs to Ns = 4 spin average (see Eq.(4)). The arrows connect the closest dots (in the sense of 
"magnetic" Euclidean distance from Eq.(5)). The columns of intra-dot structures correspond to integer multiplies 
of At = 0.01 units. 



4 




1010 1020 1050 1110 1160 1200 



1300 



Figure 6: The detailed view on reversal in the ART 
representation for sufficiently small p — 0.8 showing 
the process of formation and annihilation of intra-dot 
solitons (for the same parameters used in previous 
picture) . 

5. CONCLUSIONS 

The magnetization reversal has been simu- 
lated for superlattice model of ultra-dense mag- 
netic dot array. Results of ART classification 
demonstrated that details of reversal are deter- 
mined by the interplay of exchange and magneto- 
static couplings, inter-dot domain wall couplings 
and finite-size effects. The consequence of peri- 
odic exchange-coupling cuts is the formation of 
intra-dot soliton-antisoliton pairs. 

The neural network is a new type of model 
which allow studying of systems with many in- 
teracting pieces. We have shown that ART dia- 
grammatic viewpoint is helpful in thinking how a 
dot switches from one state to another within the 
large time intervals, and how to extract the typi- 



cal channels of intra-dot evolution. The computa- 
tional experience creates also believe that advan- 
tages of ART paradigm should be valued namely 
in 3D case, where visual classification meets hu- 
man recognition bounds. 
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